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ABSTRACT 


Recent observations of GeV/TeV photon emission from several X-ray binaries 
have sparked a renewed interest in these objects as galactic particle accelerators. 
In spite of the available multi-wavelength data, their acceleration mechanisms are 
not determined, and the nature of the accelerated particles (hadrons or leptons) 
is unknown. While much evidence favors leptonic emission, it is very likely that 
a hadronic component is also accelerated in the jets of these binary systems. The 
observation of neutrino emission would be clear evidence for the presence of a 
hadronic component in the outflow of these sources. 

In this paper we look for periodic neutrino emission from binary systems. 
Such modulation, observed in the photon flux, would be caused by the geometry 
of these systems. The results of two searches are presented that differ in the 
treatment of the spectral shape and phase of the emission. The ‘generic’ search 
allows parameters to vary freely and best fit values, in a ‘model-dependent’ search, 
predictions are used to constrain these parameters. 

We use the IceCube data taken from May 31, 2007 to April 5, 2008 with 
its 22-string configuration, and from April 5, 2008 and May 20, 2009 with its 
40-string configuration. For the generic search and the 40 string sample, we 
find that the most significant source in the catalog of 7 binary stars is Cygnus 
X-3 with a 1.8% probability after trials (2.1 o one-sided) of being produced by 
statistical fluctuations of the background. The model-dependent method tested 
a range of system geometries — the inclination and the massive star’s disk size 
- for LS I +61° 303, no significant excess was found. 



1. Introduction 


Recently, the observation of Very High Energy (VHE) 7 -ray emission from several 
high mass X-ray binaries has established a new subclass of VHE or “ 7 -ray-loud” binaries. 
Gamma-rays are particularly interesting for neutrino searches because they can both 
produced by pion decays in an hadronic scenario. Three such systems, PSR B1259-63, 

LS 5039, and LS I +61° 303, have been identified as persistent TeV 7 -ray emitters, while 
two others, Cygnus X-l and HESS J0632^057, are possible candidates. PSR B1259-63 
(Aharonian et. al. 2005a) is formed by a B2Ve star orbited by a young 48 ms pulsar (Tavani 
k Arons 1997) both exhibiting a strong wind. As observed in (Aharonian et al. 2005a), its 
VHE emission could come from Inverse Compton scattering on shock- accelerated leptons 
from the interaction zone between the pulsar and wind from the star though a hadronic 
interpretation cannot be excluded. On the other hand, the driving factor of the VHE 
emission in Cygnus X-l, most probably a black hole orbiting a super-giant 09.7 star 
(Ziolkowski 2005), could be the interaction of the black hole with the strong stellar wind of 
the star. However, there has been no other evidence for steady VHE emission of Cygnus 
X-l, though a VHE flare of about lh was observed (Albert et al. 2007). 

LS I +61° 303 remains a mystery even after four decades of observations over a wide 
range of wavelengths, from radio (Gregory et al. 1979), soft and hard X-ray (Harrison et al. 
2000; Sidoli et al. 2006), GeV (Abdo et al. 2009a) and TeV photons (Albert et al. 2009; 
Acciari et al. 2009). The best measurement of its period, P, = 26.4960 ± 0.0028 d, comes 
from radio data (Gregory et al. 1999) with the orbital zero phase taken at JD 2443366.775 
(Gregory k Taylor 1978), but the same modulation has also been detected in other 
wavelengths, notably in the GeV/TeV band emission (Abdo et al. 2009a; Albert et al. 
2009). The optical component, a massive, rapidly rotating BOVe main sequence star loses 
mass through a strong stellar wind thought to be formed by a fast, low-density polar wind 
and a slow, high-density equatorial decretion disc (Romero et al. 2007). Yet, the nature of 



the compact companion cannot be constrained to exclude either a neutron star or a black 
hole (Sierpowska-Bartosik & Torres 2009). It remains open if the VHE emission stems 
from a microquasar scenario (Bosch-Ramon et al. 2006b) or a pulsar scenario (Dubus et al. 
2010). A similar argument (Paredes et al. 2000; Casares et al. 2005; Ribo et al. 2008) can 
be made for LS 5039 which was discovered in the TeV 7 band shortly after PSR B 1259-63 
(Aharonian et al. 2005b). As a consequence of our poor knowledge of the accelerator, the 
composition of the outflow is unknown. Existing observations cannot exclude the presence 
of hadrons in the pulsar wind. The observation of multi-TeV neutrinos would unequivocally 
prove the existence of hadronic acceleration. 

The paper is organized as follows. In Section 2 we describe the IceCube observatory 
and the data taken with two detector configurations (Achterberg et al. 2006). The analysis 
method is described in Sec. 3. The underlying hypothesis of the adopted likelihood method 
is that the neutrino emission is periodically modulated due to the geometry of the considered 
X-ray binaries. Neutrinos would be produced by a beam of hadrons accelerated by the 
compact object and interacting with the matter of the massive star and its atmosphere. The 
periodic modulation would be connected to the rotation of the system. This modulation is 
observed in photons from radio to X-ray, and in the VHE. While the signal is modulated 
in time, the background of atmospheric and muon neutrinos is random in time. Hence, it 
is estimated scrambling data in time - each scrambled data set is equivalent to a sample 
of atmospheric background. In the sections that follow, the search method, its expected 
discovery potential, and results are illustrated. To avoid any bias toward discovery, each 
search has been performed in a blind fashion by defining cuts before looking at the time (or 
right ascension) of the final event sample. The post-trial p- values are robust, since we use 
data themselves to estimate the background. 

Section 4 describes the first search, called hereafter ‘generic’, since it incorporates 
only minimal assumptions regarding the neutrino spectrum: a simple power law emission 
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spectrum and a period that matches the periodicity in an observed electromagnetic 
wavelength. The phase of the neutrino flux is not constrained to match the phase of the 
observed photon emission. It is instead a free parameter of a fit to account for the fact 
that photons can be absorbed when the accelerator is behind the large star of the binary 
system while neutrino production can be enhanced if enough matter is crossed. The search 
catalog consists of 7 galactic binary stars in the Northern sky. The selected objects are 
considered as microquasars in Distefano et al. (2002), where their expected emission of 
neutrinos is calculated. That paper is not specifically on the periodic emission of these 
sources, nonetheless the objects considered there are promising neutrino emitters for which 
radio observations allow identification of jet parameters such as the Lorentz factor and the 
luminosity of the jet. All the considered sources are located in the Northern Hemisphere, 
since in this region IceCube is more sensitive to sources in the TeV-PeV region (Abbasi 
et al. 2011). 

A model-dependent search was also carried out (Section 5) for LS I +61° 303. The 
serach relies on phenomenological predictions for the spectral neutrino light curve (SNLC). 
This places constraints on the nuisance parameters of this search, e.g. the spectral index 
of the neutrino emission and its phase. For LS I +61° 303, several time-integrated flux 
predictions that would lead to a detection of neutrinos in an IceCube-like detector exist 
(Distefano et al. 2002; Christiansen et al. 2006). and some authors predict a SNLC (Torres 
& Halzen 2007; Orellana & Romero 2007). The prediction used here is based on the 
emission model first mentioned in Neronov & Ribordy (2009) and developed into an SNLC 
in Chernyakova et al. (2009). In this model, the dense stellar disk that is formed by the 
outflow of the massive Be star is assumed to be oriented perpendicularly to the observer’s 
line of sight, so that the orbital phase period in which the compact object stays behind 
the massive star is maximal. The compact object’s relativistic outflow delivers high energy 
protons that are only weakly deflected in the system’s magnetic field — literature gives 
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estimates between 0.2 — 0.35 G (Chernyakova et al. 2006) and 1G (Bosch-Ramon et al. 
2006a) — so that their Larmor radius at several PeV is comparable to the size of the system. 
The interaction of these protons with the dense stellar disk subsequently emits a neutrino 
flux into a cone that is aligned with the axis of the compact object - Be star system. 
The resulting neutrino flux can only be observed if the compact object is eclipsed by the 
dense stellar disk with respect to the observer. The neutrino emission windows predicted 
by this geometric model occur before the periastron and are therefore not in phase with 
the emission maxima observed in the GeV and TeV bands, which peak at phases around 
0.3-0. 4 (GeV band) and 0. 7-0.8 (TeV band). The y-rav emission from secondary e ± and 7r° 
decay that accompanies the neutrino emission will be strongly suppressed due to the high 
density of soft photons and matter in the dense stellar disk, delaying the onset of the GeV 
emission peak. As a consequence, Chernyakova et al. (2009) does not give a pi’ediction for 
the neutrino flux normalization within their time-dependent emission model. 

2. The IceCube data 

Construction of the IceCube Neutrino Observatory started with a first string installed 
in the austral season of 2005/2006 (Achterberg et al. 2006) and was completed in December 
2010. It is composed of an array of 86 strings with a total of 5160 Digital Optical Modules 
(DOMs) instrumented between a depth of 1.450 and 2,450 m. The deep ice detector is 
complemented at the surface by the extensive air shower array IceTop. Each DOM consists 
of a 25 cm diameter Hamamatsu photomultiplier (Abbasi et al. 2010a), electronics to digitize 
the PMT output voltage (Abbasi et al. 2009b) and all in a spherical, pressure- resistant glass 
housing. IceCube observes the Cherenkov light emitted by relativistic charged particles 
produced in high-energy neutrino interactions. The PMT signals are digitized and the 
charge and arrival time of photons measured. The data taking and performances of the 
detector are described in Abbasi et al. (2011). IceCube triggers primarily on down-going 



muons: in the 40-string configuration at a rate of about 950 Hz and in the 22-string 
configuration at about 350 Hz. During the austral summer the atmosphere above the South 
Pole gets warmer and thinner and the probability of pions generated in cosmic ray air 
showers to decay rather than interact increases (Tilav et al. 2009), causing the trigger rate 
to vary by about ±10%. At final level of analysis, dominated by upwardgoing atmospheric 
neutrinos, events come from many different directions. Hence, temperature effects average 
over a wide terrestrial region and the seasonal modulation is only a few percent. In any 
case, this variation has a period of one year, much longer than any period considered in this 
search. 

The searches presented here used two data samples. The data taken with the 22-string 
configuration have a livetime of 275.7 days, or a total livetime of 89% of the operation 
period from May 31, 2007 to April 5, 2008 (MJD 54251-54561). The sample is described 
in Abbasi et al. (2009a), and consists of 5114 events, which are mostly neutrino induced 
upward-going muons with declinations from -5° to +85°. The deadtime is mainly due to 
test and calibration runs during and after the construction season. The 22-string data 
sample has been unblinded for discovery for the generic search, and for the model-dependent 
search it was used for establishing the methodology and provide an a posteriori p-value. 
Both methods were then applied on data taken with the 40-string configuration for 
discovery. The lifetime of the 40-string data used in analysis is 375.7d which is 92% of 
the nominal operation period from April 5. 2008 to May 20, 2009 or Modified Julian Day 
(MJD) 54561-54971. The handling and processing of the data to obtain the final neutrino 
candidate event sample are fully discussed in Abbasi et al. (2011). The final 40-string 
sample contains 36.900 atmospheric neutrino and muon events distributed over the whole 
sky, of which 14,121 events are up-going, and 22,779 events are down-going. The median 
angular resolution for the final samples of the analysis for energies larger than 10 TeV 
is below 1°. About 50% of the reconstructed muon events produced by neutrinos with 
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energy between 10-100 TeV (1-10 PeV) from a direction are inside 1° (0.6°). The selected 
atmospheric neutrino events are in the TeV-PeV range in the up-going region, while in 
the down-going region a high energy sample of PeV-EeV atmospheric muons is selected. 
This difference in sensitivity is due to the need of reducing the large number of downgoing 
muon events that not only would mask the neutrino signal but that are too numerous to 
apply the time-consuming likelihood method to build the full skymaps. For this scope, a 
zenith-dependent energy proxy cut has been implemented to reduce the muon event rate to 
roughly that of the atmospheric neutrino background used for the point-source search in 
the northern hemisphere. This energy proxy uses the density of photons along the muon 
track due to stochastic energy losses of pair production, bremsstrahlung and photonuclear 
interactions which dominate over ionization losses for muons above 1 TeV. The achieved 
energy resolution is about 0.3 in log 10 of the muon energy in the detector between 10 TeV 
and 10 5 TeV. The muon neutrino flux upper limits at 90% CL for time-independent searches 
(depending on declination) are between E 2 dN/dE ~ 3 — 20 x 10 -12 TeV cm -2 s -1 in the 
northern sky where the sources considered in this paper are located. 

The 22 and 40-string data samples used in this paper were also used to look for bursting 
neutrino sources in (Abbasi et al. 2010b) where the stability of the data taking is discussed 
in detail. Azimuthal geometry effects of the 40 and 22-st.ring IceCube detectors (due to the 
fact that they are more elongated in one direction than in others) and rotation of the Earth 
interfere constructively for source periods with a multiple of one half of a day. which is not 
the case for any of the source periods tested. 

The limits in this paper were produced assuming a flux of only muon neutrinos at 
the Earth with simulated energies from 10 8 to 10 19 eV. For standard neutrino oscillations 
over astronomical distances (Athar et al. 2000). equal fluxes of all neutrino flavors at 
the Earth are expected from a source producing neutrinos via pion decay with a ratio of 
v e : Vp : v T = 1:2:0. For the assumption of equal fluxes of muon and tau neutrinos at the 
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Earth, the resulting upper limits on the sum of both fluxes are about 1.7 times higher than 
if only muon neutrinos are considered (Abbasi et al. 2011). This sets better limits than the 
expected factor of two due to oscillations if no tau neutrinos would be detectable. This 
is due to the tau decay channel into muons with a branching ratio of 17.7%. In addition 
to this, tau leptons with energy greater than some PeVs that may travel far enough to be 
reconstructed as tracks in IceCube before decaying. For an E ~~ 2 neutrino spectrum the 
contribution due to the detectable tau neutrino flux for sources at the horizon is 10% and 
rising to 15% for sources in the Northern hemisphere. 

The main systematic uncertainties on the flux upper limits come from photon 
propagation in ice, absolute DOM efficiency, and uncertainties in the Earth density profile 
and muon energy loss. For an E ~ 2 spectrum the estimated total uncertainty is about 16% 
(Abbasi et al. 2011). They are included in the upper limits calculations following the 
method of Conrad et al. (2003) with the modification described in Hill (2003). 

3. Method 

The generic and model dependent searches presented here both use a likelihood 
method. The method was extensively described in Braun et al. (2008, 2010). The generic 
search requires no binning in the used variables, while the model dependent search selected 
events within an angular window from the source and with a cut on an energy proxy. In the 
likelihood ratio method, the data are modeled as a combination of signal and background 
populations. In the absence of detector biases, the background is randomly distributed in 
time, whereas the signal will have periodic features. Hence, the p-value of the results are 
calculated by comparing to the fraction of scrambled trials in right ascension (equivalent 
to time) which have a larger likelihood ratio than what found in data. The probability 
density function of the background and of the signal contains three terms: a space term, an 
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energy term and a time dependent term. The first two are described in detail in (Abbasi 
et al. 2011). The first term characterizes the signal of a point source as clustered around its 
direction with an angular spread mainly due to the angular error of the reconstruction. The 
second is used because the signal is expected to have harder spectra than the background. 
Below we will discuss each time dependent term used in the searches presented here. The 
generic search describes the shape of the periodic emission as a Gaussian. The model 
dependent search uses a rectangular window in time. In the generic method the spectral 
index is a freely fit parameter while the model-dependent method uses assumed spectra 
from models. 

We calculate the sensitivity and median upper limit at 90% confidence level using 
the method of Feldman k Cousins (Feldman k Cousins 1998). A discovery potential is 
calculated as the flux required to achieve a p- value less than 2.87xl0~ 4 * * 7 (5 a of the upper 
tail of a one-sided Gaussian) in 50% of trials. It should be noted that the threshold 
significance to claim a discovery in IceCube is set to 5a. Fig. 1 shows the sensitivity and 
the discovery potential for the two methods, together with the corresponding values from 
the time-integrated search (Abbasi et al. 2009a) . 


4. Generic periodic emission search 

Inspired by the periodic modulation of the photon emission in the TeV and GeV 7 

bands that has been observed for LS I +61° 303 (Albert et al. 2006; Abdo et al. 2009a), 

this search is extended to other microquasars and binary systems. The modulation in the 
LS I +61° 303 emission is interpreted as an indication of the absorption of 7-rays in the 
system depending on the relative position of the observer and the accelerator. We assume 
that the neutrino emission period is the same as for photons, but in this generic search we 
leave the phase as a free parameter in our method to be sensitive to possible shifts in phase 
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between neutrinos and photons. 

The best fit values for the signal fraction, spectral index, and the peak phase and 
duration of neutrino emission are obtained by maximizing the likelihood ratio. For the signal 
hypothesis, a Gaussian emission repeating each orbit is assumed. Hence our time-dependent 
pdf is: 

Ctime ^ 

* “ V^r 

where ot is the width of the Gaussian, ^ is the phase of the event and <po is the phase of 
the peak of the emission. The fit parameters are a T and <y?o- 

Comparing to the unbinned time-integrated analysis, we find that searching for 
periodicity in neutrino emission results in a better discovery potential if the duration of the 
emission a T is less than about 20% of the total period (see Fig. 1). As the time-dependent 
search adds two additional degrees of freedom to the analysis, the discovery potential is 
roughly 10-15% worse if neutrinos are actually emitted at a steady rate or over a large 
fraction of the period. For both the case of 22 strings and 40 strings (see Fig. 1), if the 
emission has a <r T of 1/50 of the period the method requires about half as many events for 
discovery than the time-independent, search. 

4.1. Results 

The seven predefined sources, listed in Table 1, were used for this search with the 
22-string sample. There was no evidence of periodicity seen for any of the sources tested. 
The most significant result for this sample was the source SS 433, with a pre-trial p- value 
of 6%, which we expect to see from at least one of our seven tested sources in 35% of 
background-only scrambled samples. 

For the subsequent 40-string configuration, the most significant deviation from the 
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background only hypothesis is for Cygnus X-3. The pretrial p- value of this source is 0.00186, 
where an equivalent best p- value from any of the sources is found in 1.8% of scrambled 
samples, so the result is compatible with random fluctuations. The peak emission is found 
to be at phase <£> 0 = 0.8, and — 0.02. The best-fit number of source events is h s = 4.28 
and spectrum is soft at % = 3.75. In a one-degree bin centered on Cygnus X-3, 4 events are 
observed with a mean background of 1.9 events. 

Fig. 2 compares the 40-string time-integrated limits to the model predictions in 
Distefano et. al. (2002) for each source. The model predicts the neutrino flux based on the 
radiative luminosity associated to the jet from radio observations in quiescent states and 
during flares the duration of which is specified in Tab. 4 in that paper. The figure shows 
limits for both the persistent and time-dependent cases for a time window similar to the 
observed flare but not coincident to it (since IceCube was not active at the time of radio 
observations noted in the paper). For the persistent case of SS 433 the model predicts more 
than 100 events during the 40-string data taking period, which is excluded at greater than 
99% condence level. It should be considered that the authors indicate in their paper that 
the model may be biassed by the fact that the source is surrounded by a the diffuse nebula 
W50 that can affect the estimate of the radio emission used in the model for SS 433. 

The main parameters on which the neutrino flux depends in this model are: the 
fraction of jet kinetic energy converted to internal energy of electrons and magnetic field, 
r ) e ; the fraction of the jet luminosity carried by accelerated protons, t] p ; and the fraction 
of proton energy in pious / ff , which strongly depends on the maximum energy to which 
protons can be accelerated. We show as an example for the case of a 3-day burst of Cygnus 
X-3 how the parameters are constrained by our result. We assume equipartition between 
the magnetic fields and the electrons and the proton component ( i] p — r ] e ) for setting a 
constraint on / # < 0.11. If equipartition does not apply, we assume = f-K.peak as given in 
Table 2 in Distefano et al. (2002) (for Cygnus X-3 f~. pea k = 0.12) and constrain rj p to be 



Source 

Period 

(days) 

p- value 
(pretrial) 

T 0 

(MJD) 

(phase) 

G’J' 

(period” 1 ) 

Reference 

Time-Dependent UL 
(TeV“ 'em" 2 s“ 1 ) 

Time- Integrated UL 
(TeV"" I cm“ 2 s~ 1 ) 

Cygnus X-3 

0.199679 

0.00186 

53760.583997 

0.819 

0.02 

(Abclo et al. 2009b) 

3.01-10 -11 

6.64-10 12 

Cygiius X-l 

5.5929 

0.080 

41874.707 

0.031 

0.02 

(Brocksopp et al. 1998) 

4.08-10 12 

7.41-10 -12 

LS I +61° 303 

26.498 

0.23 

43366.775 

0.916 

0.02 

(Acciari et al. 2008) 

1.82-10“ 11 

9.78-10 12 

GRS 1915+105 

30.8 

0.43 

53945.7 

0.502 

0.045 

(Neil ct al. 2007) 

2.57-10 12 

3.23-10“ 12 

SS 433 

13.0821 

0.35 

50023.62 

0.779 

0.02 

(Hill wig & Gies 2008) 

3.15-10 12 

3.03-10“ 12 

XTE J1 118+480 

0.169934 

0.28 

52287.9929 

0.985 

0.132 

(McClintock et al. 2003) 

7.29-10 -12 

8.18-10 12 

GRO J0422+32 

0.21214 

0.037 

50274.4156 

0.831 

0.02 

(Webb ct al. 2000) 

1.46-10 -11 

6.89-10“ 12 


Table 1:: System name, period, pre-trial p- value, and the time of zero phase for the binary systems tested, a-r is 
the fraction of the standard deviation of the best-fit Gaussian of the period of the binary system. We also include 
the reference used for the orbital information. The last columns give the upper limits of Feldman-Cousins Feldman 
& Cousins (1998) 90% confidence intervals as the normalization on an E~ 2 spectrum flux for the time-dependent and 
integrated hypotheses for the 40 strings data. The upper limits also incorporate a 16% systematic uncertainty. 
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less than 92% of r/ e . In deriving these limits we have assumed that the Lorentz factor of the 
jet is well known from radio measurements, but in many cases there is a large uncertainty 
on this parameter. Hence, our limits to the parameters of this model may have different 
implications that we cannot disentangle: protons may not be dominant in the jet, they may 
lose smaller energies into pion decay than the values considered in Distefano et al. (2002), 
or the Lorentz factor is lower than the value indicated in Table 1 in that paper. 


5. Time-dependent Search based on neutrino flux predictions for 

LS I +61° 303 

In this section, we introduce a model-dependent search for a neutrino flux from 
LS I +61° 303. In contrast to the generic method discussed above, all nuisance parameters 
of this search are fixed according to model predictions from Neronov & Ribordy (2009); 
Chernyakova et al. (2009). The model allows to make predictions for the neutrino emission 
time window from the physical parameters of the binary system. 


5.1. Phase Window Selection 

The duration of the phase where observable neutrinos are emitted is given by the SNLC 
and depends on two parameters: the inclination i of the system’s ecliptic with respect to the 
observer, and the dimensionless radius Rj D per of its dense disk. The system’s inclination 
is unknown and depends on the mass of its components. The favored mass interval of the 
optical star is M opt = 10 — 15 M e (Hutchings & Crampton 1981). With this constraint, 
Casares et al. (2005) conclude that the compact object would be a black hole for i < 25°, 
and a neutron star for 25° < i < 60°. 

In our choice of models, we consider an inclination of i = 30°, corresponding to 
the neutron star scenario with minimal mass ( 1 .4 M . ). with three different disk sizes 
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R/D pei = 1.0, 1.5, and 2.0, and an inclination of i = 10°, corresponding to the black hole 
scenario with disk size R/D peT — 1.0. RjD pcr is equal to and larger than unity, in agreement 
with the limiting values discussed in (Sierpowska-Bartosik & Torres 2009). As discussed 
above, the time-dependent model lacks a more detailed SNLC prediction. Therefore, we 
consider a flat neutrino emission, i.e. constant within the allowed orbital phase range and, 
given the lack of a precise understanding of the acceleration mechanisms at the source, we 
consider simple power law neutrino spectra. Two spectral indices of j s = 1.5 or 2 were 
chosen: the latter represents Fermi acceleration-based models, whereas the former covers 
models with stronger acceleration mechanisms, e.g. acceleration in large scale electric fields. 
In summary, each model M = M 7s ,i,/j/D per is characterized by the three physical parameters 
i , R/D per and y s . 

The bin size T and cut on the minimum reconstructed energy E cut are optimized 
according to the model detection potential (MDP) (Hill et al. 2005) by minimizing a 
hypothetical flux for a 5 a discovery in 50% trials. Table 2 summarizes the results of the 
optimization for the 22-string data. The discovery potential is expressed in terms of the flux 
normalization <I>o for a simple power law neutrino emission at the source. For this specific 
model, this approach is more sensitive than the generic one, although the differences are 
small (around 5% for the discovery potential, and 10% for the sensitivities). 

The search was optimized for the 22-string data set in a blind fashion but treated 
as an a posteriori search. Therefore, a p- value can not be stated. We find the strongest 
fluctuation over background for the window Four event candidates remained after cuts on 
the time window, search bin radius and minimum energy proxy. They are all before or 
around the superior conjunction (at an orbital phase of 0.04). Based on this observation, 
an additional time window, “22-string best”, was introduced to be tested for both spectral 
indices alongside the already introduced eight models. It allows neutrino emission between 
phases 0.9 and 1.1, giving a width of olm = 0.2 for the emission time window 7 . 
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Model Af TsjiiH /£) per 

(y’minj ^min) 

"Popt Iog(EnergyProxy) cut ajwMbg 

(deg) 

Mlds 

4>o 

(TeV -1 cm _2 s -1 ) 

M 2 , 10 , 1.1 

(0.98, 1.20) 

2.6 

2.0 

1.2 

7 

4.2-10” 11 

M 2 , 10 , 1 . 5 /^ 2 , 30 , 2.0 

(0.93, 1.24) 

2.6 

2.6 

1.5 

7 

4.5-10 -11 

M 2 , 10 , 2,0 

(0.87, 1.24) 

2.4 

2.2 

1.6 

8 

4.8-10 -11 

M 2 , 30 , 1.5 

(1.01, 1.24) 

2.3 

2.6 

0.86 

7 

4.2-10- 11 

Mi, 5 , 10 , 1.1 

(0.98, 1.20) 

2.6 

2.8 

0.92 

5 

5.610 -12 

Ml, 5 , 10 , 1.5 /Ml, 5 , 30, 2.0 

(0.93, 1.24) 

2.6 

2.8 

1.3 

6 

6.2-10 -12 

Mi. 5 , 10 , 2.0 

(0.87, 1.24) 

2.6 

2.8 

1.6 

6 

6.4-10 -12 

Ml, 5 , 30, 1.5 

(1.01, 1.24) 

2.6 

2.8 

0.97 

5 

5.7-10- 12 


Table 2:: Discovery potential for 22-string sample: a model is characterized by the inclination 
of the system w.r.t. the line of sight i, the disk size in periastron distance units R/D peI 
and the simple power law spectral index 7 ,,: M la ^ R / Dpet . (ip m in , ^ max ) are the beginning 
and the end of the neutrino emission time window, X V and EnergyProxy cut the optimal cut 
parameters, and olmIi bg the scaled average background expectation. A model leads to a > 5<r 
CL detection in 50% of trials, if an average number of events ji\^ s is detected after final cuts. 
This corresponds to neutrino flux normalization <f > 0 of a simple power law with spectral index 
7 S at the source, i.e. d<f>/d E = <f > 0 {Ej TeV)~ 7s . 


5.2. 40-string Application and Result 

The search method presented so far was re-optimized and then unblinded on the 
40-string data set. The best pre-trial p-value of 0.35 was obtained for model M1C22 best. 7s =2 
with one event remaining in the final search window. Taking into account the trial factor 
for 10 simultaneous searches, the overall post-trial p-value is 65%, which is not significant. 
Table 3 summarizes the search parameters and the results for the time-dependent upper 
limit of the 90% CL Feldman-Cousins intervals. 
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Model M ystt . R / Dpec 

'J'opt log(EnergyProxy) C ut ot M fi bg 

(deg) 

71 $ 

p-value 

(pre-trial) 

T i me- Dep endent 
Upper Limit (UL) 
(TeV -1 cm“ 2 s -1 ) 

M2, 10, 1.1 

1.4 

2.4 

0.80 

0 

0.69 

5.9-10"' 6 * * * * * 12 

A^2, 10, 1.5/M2, 30.2.0 

1.3 

2.8 

0.84 

i 

0.43 

7.9- 10- 12 

M2, 10,2.0 

1.3 

2.6 

1.1 

i 

0.51 

7.410 -12 

M2, 30, 1.5 

1.2 

2.6 

0.58 

0 

0.55 

5.9-10” 12 

Ml. 5, 10, 1.1 

1.6 

3.0 

0.71 

0 

0.78 

9.5-10~ 13 

Ml. 5,10, 1.5/Ml. 8,30,2.0 

1.3 

3.0 

0.64 

1 

0.44 

1.2-10 -12 

Ml. 5,10,2.0 

1.5 

3.0 

1,0 

1 

0.56 

1.2-10 -12 

Mi. 5, 30, 1.5 

1.6 

3.0 

0.74 

0 

0.81 

9.110 - 13 

1C-22 best, -y s = 2 

1.4 

2.4 

0.70 

1 

0.36 

8.6-10~ 12 

IC-22 best, 7 S = 1.5 

1.5 

3.0 

0.56 

1 

0.35 

1.3-10” 12 


Table 3:: Search results with 40-string data: A model is characterized by the inclination of 
the system w.r.t. the line of sight i, the disk size in periastron distance units R/D per and 
the simple power law spectral index j s : M ls , lR / Dper . T and log(EnergyProxy) cut are the 
optimal cut parameters, and o>//i b g the scaled average background expectation. n s is the 
number of source events found in the respective time window after selection cuts that give the 
corresponding pre-trial p-value. The last column gives the upper limits of Feldman-Cousins 
90% confidence intervals as the normalization on an E” 2 spectrum flux. 


6. Conclusions 

The exploration of the GeV and TeV photon sky with the instruments on board the 

Fermi spacecraft and the ground-based Cherenkov telescopes has heralded the golden age 

of y-ray astronomy. The connection to neutrino astronomy is clear: high energy processes 

which cause the observed VHE emission can be responsible for the observed high energy 

cosmic rays. This implies hadronic acceleration mechanisms in astrophysical sources which 

can result in an observable neutrino flux with giant neutrino telescopes like IceCube. 

The available photon observations have made it possible to enhance the sensitivity of 
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searches for neutrino fluxes by incorporating assumptions derived from the 7 -rav data. One 
crucial development has been the formulation of time-dependent neutrino flux searches, 
postulating a connection between the time modulation of the high energy emission and the 
possible neutrino flux. This assumption has increased the sensitivity of these searches in 
comparison to their time-averaged counterparts. 

This paper presented two searches for neutrinos from objects with periodic 
photon/broadband emission. Seven X-ray binaries in the Northern Hemisphere were 
selected and tested with the generic method which scans the whole allowed neutrino 
spectrum parameter space. The most significant source in the catalog is Cygnus X-3 
with a 1.8% probability after trials ( 2 . 1 <r one-sided). The best fit of the phase is 0.82, 
close to the superior conjunction of the system at phase = 0 and close to the peak of 
Fermi-LAT observations. Comparing the time-integrated limits for each source to model 
predictions from Distefano et al. ( 2002 ), the results exclude the prediction for SS 433 at 
greater than 99% confidence level. We show that our limits can constrain the fraction of 
jet luminosity which is converted into pions and the ratio of jet energy into relativistic 
leptons versus relativistic hadrons, under some assumptions. For instance, for Cygnus X-3 
and equipartition between electrons and protons, the fraction of proton energy in pions is 
limited to about 11 %. 

The second method searched for a neutrino flux from LS I +61° 303 based on various 
predictions of the system parameters, i.e. size of the massive star’s disk and the system’s 
inclination. The most significant model prediction yielded a 35% (36%) probability before 
trials for an emission between the phases 0.9 and 1.1 with a spectral index of j s = 1.5 (2). 
The 90% CL Feldman-Cousin time-dependent, upper limit is 1.3 • 10 ~ 12 TeV -1 cm -2 s -1 
(8.6-10 -12 TeV -1 cm -2 s -1 ). Due to the geometric nature of the underlying model, no 
statement can be made with respect to the validity of the model predictions. However, 
these upper limits can guide further developments of models. LS 1 +61° 303 was also 
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included in the catalog of the generic search, with a pre-trial p-value of 23%. Its best fit 
phase of <p 0 = 0.92 agrees well with the time window that gave the best p-value in the 
model-dependent search, (0.9, 1.1). All of the results in this paper are compatible with a 
fluctuation of the background. 
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Fig. 1. — : Discovery potentials (5a in 50% trials, solid and dashed upper lines in the plot) 
and the average sensitivity at 90% CL calculated with the prescription in (Feldman & Cousins 
1998) (solid and dashed lower lines) for the generic search applied to the 22-string data (left) 
and to the 40-string data (right). The y-axis shows the flux normalization < F 0 for a E ~ 2 
simple power law spectrum emission from LS I +61° 303 , i.e. d$/dE = <f> 0 (£’/TeV)"" 2 . The 
x-axis shows the width a^ of the Gaussian emission normalized to its period Pi — 26.4960 d. 
Shown are the results for the time-integrated searches (Abbasi et al. 2009a, 2011) (dashed 
line) and for the generic time-dependent search (solid line). Symbols are the results for the 
model-dependent search (for comparison purposes in this plot, the rectangular time emission 
window otherwise used by the model-dependent search were converted into a Gaussian with 
the same normalization). 
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Fig. 2. — : The time-integrated upper limit (UL) at 90% CL on the number of events is 
compared to the expected number of events for model predictions according to Distefano 
et al. (2002) for specific sources for the 40 string configuration. The neutrino energy range 
used to calculate the total number of events is 10 n — 5 x 10 14 eV, comparible to what was 
assumed in the model. For non persistent but flaring sources, the parameters of the model 
were estimated for flares observed before IceCube construction. Hence the time-dependent, 
upper limits are calculated averaging over a duration equal to the model flare during 40- 
string data taking (indicated as MUL in the legend). LS I +61° 303 is assumed to be a 
periodic flaring source in a high state during 26% of the orbit. 
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